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Abstract. A large fraction of the baryons today are pre- 
dicted to be in hot, filamentary gas, which has yet to be 
detected. In this paper, we use numerical simulations of 
dark matter and gas to determine if these filaments and 
groups of galaxy will be observable be XMM. The simu- 
lated maps include free- free and line emission, galactic ab- 
sorption, the XMM response function, photon noise, the 
extragalactic source distribution, and vignetting. A num- 
ber of cosmological models are examined as well as a range 
of very simple prescriptions to account for the effect of su- 
pernovae feedback (preheating) . We show that XMM has a 
good chance of observing emission from strong filaments at 
z ~ 0.1. This becomes much more difficult by z ~ 0.5. The 
primary difficulties lie in detecting such a large, diffuse ob- 
ject and in selecting an appropriate field. We also describe 
the range of group sizes that XMM should be sensitive to 
(both for detection and spectral analysis) , although this is 
more dependent on the unknown nature of the feedback. 
Such observations will greatly improve our understanding 
of feedback and should also provide stronger cosmological 
constraints. 

Key words: Radiation mechanisms: thermal - Methods: 
N-body simulations - Galaxies: clusters: general - Cosmol- 
ogy: large scale structure of the universe - X-ray: general 



1. Introduction 

The filamentary large scale distribution of matter in the 
universe is a consequence of the gravitational instability of 
a medium that was relatively smooth in early times, with 
the non-linear patterns we see today having evolved from 
small density fluctuations. There have been two compet- 
ing stories for the formation of structures: the (Russian) 
pancaking picture (Zel'dovich 1970) and the (Western) 



hierarchical clusteri ng picture (Peebles 1980J ). The "cos- 
mic web" paradigm (Bond et al 1996) reconciles these two 
scenarios, showing that the final-state is actually present 
in embryonic form in the overdensity pattern of the ini- 
tial fluctuations, with nonlinear dynamics just sharpening 
the image. This has several immediate applications to ob- 
servations, in particular: superclusters are predominantly 
cluster-cluster bridges and the strongest filaments exist 
between aligned clusters. 

The study of galaxy clusters has undergone spectacular 
developments in the last decade, ranging from the galaxy 
content to the detailed physics of the hot intra-cluster gas 
(ICM) and also addressing the question of the large scale 
distribution of clusters. Conversely, intermediate density 
structures, that we would define as density contrasts of the 
order of Sp/p ~ 4, i.e. loose galaxy groups and filaments, 
have captured less attention so far. The main reason being 
their comparatively humble appearance in the optical or 
X-ray sky. 

However, both theory and simulations indicate that 
a better knowledge of the filaments will improve our 
understanding of cosmology. Many current models sug- 
gest that a large fraction of the baryons at low redshift 



should be in filaments (e.g. Miralda-Escude ct al. 1996 
Cen fc Ostrikcr 1999] ). Heated to a temperature intcrme- 
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diate between that of hot clusters and cool voids, this 
would imply that presently most of the baryons remain 
undetected. Moreover, the overall cosmic network - best 
underlined by the numerous groups and filaments - may 
show different topologies depending on the value of the 
cosmological parameters and the nature of dark matter. 

Filaments and groups are superb laboratories for 
studying the effects of galaxies on their environment and, 
in turn, the impact of this environment on galaxy evolu- 
tion. The velocity dispersion of groups and filaments (~ 
400 km/s) is larger than the internal velocities of galax- 
ies themselves, thus interactions are frequent and the in- 
tragroup gas should have very similar properties to the 
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neighbouring filamentary gas. Feedback from galaxy/star 
formation is expected to play a major role - - much 
more, probably, than for massive clusters where gravi- 
tational physics may be sufficient to explain the bulk of 
their properties. Also, many quasars are thought to be 
in groups and may have significant effects on their en- 
vironment. In the various observed X-ray/optical correla- 
tions, groups tend to show a larger dispersion than clusters 
(e.g. Mulchaey fc Zabludoff 1998|) . There are also some 
systematic difference in the general properties of clusters 
and groups ( Rcnzini 1997 ). These results certainly suggest 
that galaxies play an important role in groups and fila- 
ments, but this process is very poorly understood. This is 
one of the compelling arguments in favor of deep and ex- 
tended X-ray surveys: a large sample of groups (free of pro- 
jection effects or of any optical selection biases, which can 
be very stringent e.g. "compact" groups) would enable the 
statistical study of their temperatures and abundances to 
be further compared with their optical properties and the- 
oretical predictions. Since groups occupy the intermediate 
range in the cosmic mass function — between clusters and 
galaxies — a better understanding of their global proper- 
ties would provide an alternative view of the missing link 
between quasi-linear and extremely non-linear systems. 

Our observational knowledge of the truly filamen- 
tary medium is even slimmer than for groups, border- 
ing on non-existent. Attempts to detect the "superclus- 
ter" or filament gas directly from its bremsstrahlung 
emission have led to very marginal or negative 
results so far (Briel & Henry 1995[ |Wang et al. 1997; 



Kull fc Bohringer 1999| ; see also |Markevitch 1999| ). Al 
though the temperature of this gas is expected to be 
around 1 keV, its low density considerably limits its poten- 
tial emissivity in the X-ray or EUV band. It is thus primar- 
ily a technical challenge to observe this warm component 
which is predicted to account for so me 20-40% of the tota l 
baryonic mass by redshift z = ( |Cen fc Ostrikcr 1999 ). 
In addition to a direct window into the physics of this 
medium, closely connected with the cooler components 
such as Lya clouds and galaxies, quantitative observations 
would put an upper limit to the contribution of filaments 
to the auto-correlation function of the X-ray background. 
Thanks to its unrivaled sensitivity, XMMtl is to open a 
new era in the [0.1-10] keV range. One aspect of the mis- 
sion is focussed on high resolution spectroscopy, the other 
main one being the observation of faint extended diffuse 
emission. In the EPIC imaging mode at 1 keV, XMM has 
a Half Power Diameter of 13" on-axis, a collecting area 
of 2600 cm 2 and an energy resolution of 20%. Beside the 
detailed study of massive clusters of galaxies, XMM is 
thus ideally suited to search for low surface brightness ex- 
tended structures such as very distant galaxy clusters (at 
redshifts substantially greater than one, if they exist!) or 
groups and filaments out to a redshift close to unity. 



XMM instruments are now fully calibrated and the 
launch is scheduled for the end of the millennium. It is 
thus timely to investigate carefully which of the issues re- 
lated to intermediate density structures can be efficiently 
addressed by XMM. This will have a direct impact on 
the preparation of dedicated observations as well on their 
analysis. In order to provide quantitative arguments, we 
have simulated XMM images. The X-ray emissivity of a 
filamentary region at a range of redshifts was computed 
from high resolution hydrodynamical simulations com- 
bined with a detailed plasma code and folded with the 
XMM spectral and imaging responses. We present here the 
XMM simulated images and further discuss the physics 
that could be extracted from them. 

The paper is organized as follows. In the next sec- 
tion, we present the hydrodynamical code and the various 
cosmological models and epochs that were chosen for the 
the simulations. Section 3 describes how the X-ray images 
were produced. The general discussion on the detectabil- 
ity of groups and filaments is given in Section 4 as well as 
practical hints for future large surveys of groups. In the 
conclusion we present further observational prospects in 
connection with other wavelengths. 

2. The AMR simulations 

2.1. Basic ingredients of the simulations 

In order to generate realistic XMM images, we have per- 
formed a series of cosmological simulations incorporating 
both dark matter and baryonic gas. The method used 
is an adaptive mesh refinement (AMR) algorithm that 
models the dark matter with particles and the gas with 
a series of meshes. More compl ete descriptions of the 
method are available e lsewhere ( Bryan fc Norman 1997 ; 
Norman fc Bryan 199£ ), but we briefly outline it here. 
An initial grid was set up at high redshift with density 
and velocity fields drawn from a power spectrum appro- 
priate for the cosmology of interest. The particles were 
initialized using the Zel'dovich approximation. The equa- 
tions of hydrodynamics were solved on the mesh using the 



picccwisc parabolic method (Colella & Woodward 1984; 
Bryan et al. 1995|). Since the mesh is fixed, gravitational 
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collapse inevitably causes structures to shrink to sizes 
comparable to the mesh spacing, causing inaccuracies. 
AMR solves this problem by identifying regions which 
need enhanced resolution and laying down a finer mesh 
over those volumes. The solution is interpolated onto this 
finer grid — which typically has cell sizes half as large — 
and an improved solution is computed there. The process 
can be repeated as necessary and a hierarchy of grids is 
generated, with cell sizes varying as small as necessary 
(up to a predefined limit in order to limit the cpu time 
required). 

In Table |], we list the parameters of the simulations 
performed. They are drawn from three different cosmolo- 
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gies: 1) a standard cold dark matter (SCDM) model, in 
order to link with previous work on this model, 2) a cos- 
mological constant-dominated universe (ACDM), which 
agrees with the majority of current observational data, 
and 3) a tilted-CDM (tCDM) model. Together, these three 
choices should provide a good range in the currently viable 
parameter space. 

In addition to changing the cosmology, we also ex- 
perimented with additional physical processes in the 
simulations. In some simulations we included a sim- 



ple equilibrium cooling-curve ( Sarazin fc White 1987 ) (al- 
though we should note that the spatial resolution 
of these simulations, 50 hT 1 kpc, were not sufficient 
to fully resolve the cooling instability). The physical 
state of the intergalactic medium is intimately linked 
with galaxy formation (through galactic outflows) and 
hence star formation and supernovae. In fact, simu- 
lations of clusters (e.g. Navarro, Frenk fc White 199GJ ; 



Bryan & Norman 1998) show that gravity and adiabatic 
gas physics alone do not correctly predict the luminosity- 
temperature relationship. Instead of the observed relation 
which is approximately Lx ~ T 3 , they predict Lx ~ 
T 2 . The difference is most likely due to feedback from 



stars (Kaiser 1991; Balogh, Babul fc Patton 1999), which 



those simulations did not include. We would like to take 
into account the effect of feedback here, unfortunately 
the input of energy from supernovae and stellar winds 
i s technically difficult to include properly in simulations 
( Metzler fc Evrard 1994 ). In place of a full treatment, we 
attempt to include this effect by pre-heating the gas, effec- 
tively setting a minimum entropy (see also Navarro, White 
and Frenk 1996). This does not have a large effect on the 
hot gas in the outskirts of clusters which, due to shocking, 
naturally has a high entropy. However, it does produce a 
larger core radius, decreasing the central gas density com- 
pared to a simulation without pre-heating. This decreases 
the luminosity of small clusters and groups (since the X- 
ray emissivity is proportional to the square of the density), 
and can reproduce the observed luminosity-temperature 
relationship. Since it is not clear what this initial entropy 
should be, we tried two cases: a mild feedback in which 
the initial gas temperature is set to 4 x 10 7 K at z = 20 
(and would adiabatically cool to 9 x 10 4 K by z = for a 
fluid element which follows the cosmic mean density with- 
out shocking), and a more extreme version in which the 
initial temperature is 5 times higher. The impact on the 
Lx — T relationship of introducing galaxy feedback in the 
simulations is shown in Fig. pi. This form of "feedback" is 
successful in reproducing the observed relation. Clearly, a 
simple pre-heating of the gas does not produce the same 
energy distribution as the feedback from millions of mas- 
sive stars; however the extremes adopted here (no feed- 
back at all contrasted with heating the gas everywhere) 
will at least indicate which diagnostics are sensitive to the 
feedback and by how much. 
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Fig. 1. The introduction of "mild" feedback (see text) in 
the simulations results in a decrease of the luminosity of 
small groups and, consequently, steepens the Lx — T rela- 
tionship in agreement with the observations: observations give 
~ Lx oc T 3 whereas Lx oc T 2 is expected in the case of a 
pure gravitational collapse. The observed correlation appears 
to be very tight when "cooling flow" clusters are excluded 
( Arnaud fc Evrard 1999) and shows n o obvious evolution out 
to z ~ 0.4 ( |Mushotzky fc Scharf 1997J ) 



The cosmological constant-dominated model we have 
selected is in good agreement with most available observa- 
tional data, and in particular with the shape and evolution 
of the luminosity function. The flat, tCDM model is cer- 
tainly less-favoured by the present data. It is consistent 
with the number-density of clusters and COBE measure- 
ments of the microwave background; however, it does not 
match with the many indications of a low-density universe, 
nor with the recent supernovae data implying the presence 
of a cosmological constant. We include it to indicate the 
range of variation that cosmology will introduce in our 
filament properties. The X-ray properties of the models 
are particularly important for this study and so in Fig- 
ure g, we plot the luminosity function for each model, as 
predicted from a simple Press-Schechter prescription com- 
bined with a relation between mass and temperature, cou- 
pled with the observed luminosity-temperature relation 
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(see [Bryan fc Norman 1998 for more details on this proce- 
dure). While both models agree at z = 0, observations are 
now becoming available at much higher redshifts. For ex- 
ample, the lack of evolution out to z ~ 0.8 for moderate to 



high-luminosity clusters (Rosati et al. 1998) places tighter 
constraints on the models. In Figure |2| we also plot the 
predicted luminosity function at z = 0.8. For the ACDM 
model, the lack of evolution is naturally predicted; how- 
ever for the tCDM model to agree, we have to stretch the 
bounds of the observed evolution (or rather lack of evolu- 
tion) of the luminosity-temperature relation, which is the 
most uncertain part of the modeling procedure. For exam- 
ple, using a very recent determination of the evolution of 
the Lx — T relation ( Reichart, Castander fc Nichol 1999 ), 
the tCDM model disagrees at the 2 a level. Note, how- 
ever, that the observational luminosity function at high 
redshift is still somewhat uncertain (see for example 



Reichart et al. 1999, where a deficit of very high luminos- 



ity clusters is observed). 
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Fig. 2. The predicted luminosity function of the ACDM (top) 
and tCDM (bottom) models at z = and z — 0.8. In the 
bottom panel, we also show the effect of varying the amount 
of redshift evolution in the observed Lx — T relation by one 
and two sigma. Since observations seem to indicate little or 
no evolution in the luminosity function, the tCDM model is 
marginally inconsistent. 



2.2. The extracted filament 

Since we are primarily interested in the filament environ- 
ment, we specifically selected such a region. This was done 
by initially simulating — at low resolution — a cube 100 
/i _1 Mpc on a side. The output was examined, and a vol- 
ume 18 x 23 x 33 h~ 3 Mpc 3 was selected. This region lay 
between two large clusters and was the strongest filament 
in the box. In order to better sample the initial condi- 
tions, the simulation was rerun with an addition mesh of 
size 28x34x50 placed over the region. As many particles 
as grid points were used. The adaptive meshing was then 
turned on and a relatively high resolution simulation of 
this region was performed. The rest of the 100 h^ 1 Mpc 
cube was again simulated at low resolution, which is suf- 
ficient to provide the necessary gravitational tidal field. 



3. Simulations of the XMM images 

From the temperature and density cubes, we have calcu- 
lated for each simulation cell, the X-ray plasma emissivity 
using a redshifted Raymond-Smith code having a metal- 
licity of 0.3 solar. The choice of this latter value - which is 
comparable to that of galaxy clusters, although the total 
production of heavy elements by group galaxies is cer- 
tainly lower - has been guided by the fact that mixing 
and enrichment are probably more efficient than in clus- 
ters, the medium being much less dense. Then, assuming 
a given exposure time, we computed the photon numbers 
to be received from each cell by XMM taking into account 
the complete instrumental response (a composite matrix 
made of the responses of the 2 MOS CCDs and of the PN 
CCD in imaging mode folded with the X-ray telescope re- 
sponse). This was done for the [0.4-4] energy band, best 
suited for detecting structures having temperatures of the 
order of 1-2 keV; we do not consider the photons below 
0.4 keV as diffuse emission is too much affected by galac- 
tic absorption. In the calculation, we assumed a galactic 
hydrogen column density of Nu = 5 10 20 cm 2 . (Note that 
assuming a metallicity of 0.1 instead of 0.3 would decrease 
the count rate of a 2 keV group or of 1 keV filament by 
10% or 40% respectively in the [0.4-4] keV band). 

For each set of AMR simulations we thus obtained 
a "photon cube" which was in turn projected along the 
line of sight, i.e. the x axis. The resulting photon image 
has a pixel size equal to the apparent size of the AMR 
cell seen at the simulation redshift (Table 2, step la); 
Poisson photon noise was then added to each pixel (step 
lb). This is typically what an XMM wide angle survey 
would look like once the pointings have been placed side 
by side, all instrumental effects corrected, extragalactic 
point sources removed and the resulting image rebinncd 
to increase the signal (note however, this is an ideal case 
without diffuse background). In a second step, we have 
simulated in more detail some individual pointings, fo- 
cussing on a given area of the projected photon cube as 
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Table 1. AMR simulation models - comoving cell size: (50 lcpc/h) 3 



Model 


0"8 


h 


^ ^matter 


«A 


n b 


n 


cooling? 


"feedback"? 


CDM 


0.95 


0.6 


0.4 


0.6 


0.06 


1.0 


no 


no 


ACDMfl 


0.95 


0.6 


0.4 


0.6 


0.06 


1.0 


no 


mild 


ACDMc 


0.95 


0.6 


0.4 


0.6 


0.06 


1.0 


yes 


no 


ACDMcfl 


0.95 


0.6 


0.4 


0.6 


0.06 


1.0 


yes 


mild 


ACDMf2 


0.95 


0.6 


0.4 


0.6 


0.06 


1.0 


no 


extreme 


tCDM 


0.6 


0.5 


1.0 


0.0 


0.06 


0.81 


no 


no 


SCDM 


0.6 


0.5 


1.0 


0.0 


0.06 


1.0 


no 


no 



Table 2. Steps in the simulation of the XMM images 



Step 


Operation 


Oa 
Ob 


AMR cube + plasma code — > X-ray emissivity cube 
Oa + XMM response + Exp. time — > photon cube 


la 
lb 


projected photon cube 
lb + Poisson noise 


2a 
2b 


a single XMM field of view from la 
2a + extragalactic point sources 



would be seen by the circular XMM field of view of ~ 30 
arcmin. Practically, this consists of extracting a given sky 
region from step la, rebinning it with a pixel size of 4 arc- 
sec and then applying the vignetting function allowing for 
Poisson statistics in each pixel (step 2a) . Finally, we have 
introduced the extragalactic LogN-LogS source distribu- 



tion (Hasinger et al 1998) assuming that all sources have 



a power law spectrum with a photon index of 2. Images of 
these pointlike sources have been convolved by the mean 
XMM/EPIC point spread function as a function of off- 
axis angle and randomly added to the filament sub-image 
together with a constant diffuse background rate of 10 -5 
count/s/pixel(4") (step 2b). The instrumental character- 
istics used in the simulations are from the ground-based 
calibration phase (e.g. Dahlcm fc Shartel 1999 ). 



Table 3. Simulated image parameters 



step 


redshift 


pixel size 
kpc & arcsec 


image size 
pixels & arcmin 


1 
1 


0.1 
0.5 


25/h 19" 
100/h 19" 


952 x 1392 295' x 432' 
238 x 348 74' x 108' 


2 
2 


0.1 
0.5 


~ 5/h 4" 
~ 20/h 4" 


512 x 512 = 34' 
512 x 512 = 34' 



In order to directly compare the filament emissivities at z = 0.1 
and z = 0.5, we have rebinned the z — 0.1 image obtained in 
step la with a pixel size 4 times smaller than the original simu- 
lation cell of 100/h kpc. In the calculations, we have neglected 
the ~ 10% difference in angular distance between the open and 
closed models at z — 0.5. 



4. Results 

In order to study the effects of cosmology, distance and 
input physics on the detectability by XMM of groups and 
filaments, we display the simulated images of the entire 
filament, prior to including the back/foreground popula- 
tion of X-ray point sources. We first discuss the visibility 
of the z = 0.5 case, then the z = 0.1 one, as the larger ap- 
parent size of the diffuse structures is likely to make their 
detection more problematic, despite their proximity. 

4-1. Global images at z = 0.5 

Fig. 8 shows the predicted XMM images for our three cos- 
mological models, namely: standard CDM, tilted CDM 
and ACDM. As expected, the difference between open 
and closed models is striking. Conversely, the tCDM and 
SCDM models appear to be very similar in terms of X- 
ray emissivity, at least for the group population. This 
seems quite reasonable, since both were normalized at the 
8 Mpc/h scale. In Fig. 9, we display four realizations of the 
ACDM model calculated under various heating or cooling 
assumptions as given in Table 1, to be compared to the 
ACDM model of Fig. 8. 

Considering first the effect of cooling alone (ACDM & 
ACDMc models), groups appear substantially smaller in 
spatial extent when cooling is at work. This is because the 
cooling has reduced their pressure support, which keeps 
them extended. Moreover, the ACDM & ACDMc den- 
sity maps shown in Figs. |] & [| indicate that the model 
with cooling has more small clumps, but most of them do 
not show up in the X-ray images: these are galaxy-sized 
objects, unresolved by the simulation. In these clumps, the 
gas can cool down to 10 4 K, which prevents most of them 
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from cmmiting in the X-ray band. In the larger groups 
however, the gas in the center is still hot enough to ra- 
diate and much denser than in the no-cooling case which 
makes their X-ray profile significantly more peaked. 



Fig. 3. Density cube for the z — 0.5 ACDM model shown for 
an isosurface of 4 10 -4 atom cm -3 



4-2. Individual pointings at z = 0.5 

Figures 10 & 11 show the two individual pointings indi- 
cated on Fig. 9. Although groups as small as 10 14 M Q (or 
smaller depending on the model) will be visible for inte- 
gration times as short as 10 ks, the diffuse filamentary gas 
seems not to be detectable even with 200 ks exposures. 
This is essentially due to the very low density of the fil- 
ament (Fig. ||) and the high background level caused by 
source confusion. It is interesting to note that in all models 
the intra-filament gas is expected to have a temperature 
around 1 keV as can be appreciated on Fig. o 



Fig. 5. Density cube for the z — 0.5 ACDMcfl nrodel, shown 



Fig. 4. Density cube for the z — 0.5 ACDMc model shown for 
an isosurface of 4 10 -4 atom cm -3 



As discussed in Sec. 2.1, preheating raises the gas en- 
tropy prior to gravitational collapse, preventing high cen- 
tral densities by the formation of a larger core. This is 
best seen when considering the "extreme" feed back model 
(ACDMf2 Fig. 9): clusters are much more extended, and 
small objects are virtually absent. Group central densities 
are significantly lower than in the ACDM models, with 
only three clumps reaching the limiting density of 4 10 
atom cm~ 3 (to compare with Figs. || & [|). Some of the 
gas which did not succeed in collapsing is now trapped in 
the filamentary structures which explains the higher emis- 
sivity of the diffuse region compared to any other models. 
Paradoxically, the average temperature in the filament is 
somewhat lower in ACDMf2 than in ACDM; this is due 
to the absence of shock-heating in the pre-heated filament 
gas. These effects are less pronounced if mild feedback is 
assumed. 

Finally, we considered a mixed case, combining the ef- 
fects of feedback plus cooling (ACDMcfl), likely to be the 
most realistic. This shows very little difference from that 
with feedback alone (contrary to the striking discrepancy 
between ACDM and ACDMc) suggesting that here the 
gas stays too hot and diffuse to cool. While none of these 
models captures the detailed dynamics of real stellar and 
galactic feedback, they do give an indication of the direc- 
tion and amplitude of their effects. 



Fig. 6. Temperature cube for the z = 0.5 ACDMcfl model, 
shown for an isosurface of 1 keV 



4-3. Global images at z = 0.1 

The apparent size of the same filament at z = 0.1 is very 
large (~ 5° x 7°) and it would require more than 300 point- 
ings to be entirely covered by XMM (assuming a spacing 
of 20' per pointing to minimize the loss of sensitivity due 
to vignetting). Nevertheless, as for the z = 0.5 models, we 
have calculated the number of photons to be collected by 
XMM from the entire structure for individual integration 
times of 200 ks. The results are summarized in Fig. 12 
for the SCDM and ACDM models with or without cool- 
ing/heating. Compared with the corresponding z = 0.5 
images, one sees clearly the evolution of structure in terms 
of merging: small groups are flowing along the filaments 
to reach bigger entities located at the nodes. 

4-4- Individual pointings at z = 0.1 

We have performed detailed simulations of two XMM 
fields specifically centered on filamentary segments (cf Fig. 
12). The images are presented in Fig. 13 and suggest that, 
at this redshift, very faint details of the network will be 
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detectable. These may include collapsing regions as well 
as pure diffuse filamentary emission. 

5. Prospects for XMM observations 

5.1. The u nb (photons) - Mass (groups)" relation 

In order to get a better idea of the observability of the 
groups seen in these simulated XMM images, we show 
in Fig. [7j the relation between the mass of a group and 
the number of photons observed for the ACDM simula- 
tion at z = 0.5. The number of photons measured comes 
from extracting each region above a given flux threshold 
from the photon cube. The mass is calculated by summing 
both the gas and the dark matter within a sphere whose 
mean density is about 300 times the mean density of the 
universe (see Bryan fc Norman 1998 for a more complete 
definition). A power law fit produces a slope of 0.59. Sim- 
ple adiabatic scaling laws would indicate a relation like 
M ~ L^ 75 j while our feedback simulations produce some- 
thing like M ~ L° x 5 ~ 06 . Note, however, that we do not 
directly measure the emissivity, but count individual pho- 
tons to be detected by XMM in the [0.4-4] keV band (as- 
suming Nh = 5 10 20 cm~ 2 ). The scatter in the relation is 
about 20%. 



10 14 M will produce some 100 photons in 10 ks thus, en- 
abling simple source extent characterization as can be ap- 
preciated in Fig. 11. More generally, with 10 ks exposures, 
the 3cr detection threshold for extended sources having a 
characteristic size of about 1 arcmin ( i.e. a 2x 250 kpc 
core radius seen at z ~ 0.5 — 1) is around 65 photons 
on-axis, which is 2 time s higher than for pointlike sources 
( IValtchanov et al 2000| ). 

In the previous section, we demonstrated that both 
preheating and cooling affects the detectability of small 
groups by making them either more centrally con- 
centrated (radiative cooling) or more extended (pre- 
heating). Of course it should be kept in mind that 
preheating is simply our attempt to model stellar 
feedback. More sophisticated methods are being de- 
veloped which attempt to account for the multi- 
phase nature of the flow (Wu, Fabian & Nulsen 1998; 



Teyssier, Chieze & Alimi 1998) as well as 



more com- 
plete models of galaxy formation and feedback (see for 



example Yepes et al. 1997; Huffman fc Pharasyn 1999 



Martin 1999 ). In developing and testing these methods, 
observational comparisons will play an important role. 
Therefore, a large shallow XMM survey would provide a 
unique data base to study, in a statistically complete way, 
the physics of groups. 
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Fig. 7. The relation between group mass at z = 0.5 and the 
number of photons detected in the 0.4-4 keV range from a 
simulated 100 ks XMM image (ACDM model). 



5.2. Shallow pointings 

With exposure times as short as 10 ks, XMM will reach a 
sensitivity of ~ 5 10~ 15 erg/s/cm 2 in the [0.5-2] keV band 
for pointlikc sources. Our simulations (Fig. 11) show that 
at this limiting flux, objects having masses down to 10 14 
Mq will be readily detectable out to z ~ 0.5. This in- 
cludes all clusters and much of the group population. As 
indicated in Fig. M, groups having masses of the order of 



5.3. Deep pointings 

Assuming 200 ks exposures, XMM can reach a sensitivity 
of ~ 10~ 16 erg/s/cm 2 but will then face the generic prob- 
lem of source confusion. The limiting flux where this is 
expected to occur corresponds to exposure times between 
30-50 ks. The exact level will remain unknown until the 
first XMM or Chandra in-flight observations. This is due 
to the lack of previous deep high resolution observations 
above 2 keV (ASCA had a PSF of - 3 arcmin). It is thus 
obvious that the main hurdle for studying diffuse emission 
(clusters or filaments) will be contamination by the point 
source population, rather than the limited number of pho- 
tons. Note that, to the contrary, source confusion should 
barely be a problem for Chandra (PSF ~ 1 arcsec) but 
its much lower effective area makes the task of detecting 
filaments much more difficult. 

An additional difficulty comes from the fact that the 
simulations have a finite thickness, which is not the case 
of the "observable" universe; thus, true confusion may be 
even worse than presented here because of the presence of 
fore/background filaments. 

5.3.1. Groups 

Assuming a 50% efficiency for the X-ray detection algo- 
rithm in measuring extended source fluxes and that a net 
minimum photon number of 1000 is necessary for spec- 
tral studies at low energies, Fig. [7j indicates that it will be 
possible to obtain temperature and abundance estimates 



s 
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for objects having masses down to ~ 1O 14 M0 at z = 0.5 
(i.e. some 2000 photons collected during an exposure time 
of 200 ks). It should also be kept in mind that we are 
assuming that the X-ray emission from AGN and other 
galactic sources can be identified and removed from the 
group X-ray emission. 

Such observations will be particularly important as it 
will provide new insights into the physics of low mass clus- 
ters and groups. For example, since supernovae are respon- 
sible for both adding energy and metals to the intergroup 
medium, a measurement of the metallicity helps to con- 
strain the feedback history (which is even more important 
in groups than in clusters). Determinations of the baryon 
fraction can also cast light on the strength of feedback and 
the nature of group evolution. 

5.3.2. Filaments 

The detection of filamentary emission (and not simply of 
"supercluster gas" ) is a real challenge for the next genera- 
tion of X-ray satellites. With its unrivaled sensitivity and 
large field of view, XMM is a priori the best suited instru- 
ment. Our detailed simulations indicate that this will be 
very difficult at z ~ 0.5 whatever the assumed cosmology 
or heating/cooling processes, but possible at z ~ 0.1 for 
an open model associated with mild feedback (ACDMfl, 
Sec. 4.4). 

More generally, in this work, we have tried to gauge the 
effect of uncertainties in our understanding of the physics 
of filament (i.e. cooling and feedback). We have shown that 
cooling alone tends to decrease somewhat the luminosity 
of filaments. This occurs because gas cools into many small 
clumps, reaching temperatures which are too low to radi- 
ate in the X-ray band. However, this also produces groups 
and clusters which are much brighter than observations 
indicate (as reflected in the mismatched Lx — T relation) . 

On the other hand, in the simulation with the extreme 
form of preheating (ACDMf2), small groups are prevented 
from forming at all (in a more realistic simulation, these 
groups would form but then expel nearly all of their gas 
due to supernovae in galaxies). This causes more gas to 
end up in the filaments, improving their visibility. 

The case of mild feedback, which satisfies all the ob- 
servational constraints is intermediate in terms of fila- 
ment emission. Feedback and cooling (ACDMcfl, the most 
likely model) appears very similar to the model with feed- 
back alone (ACDMfl) in terms of emissivity. Thus, our 
prediction that filamentary emission will be detectable at 
z ~ 0.1 is likely to be robust. Note also that the varia- 
tion due to changes in input physics is much smaller for 
filaments than for groups and filaments, which are more 
sensitive to the nature and strength of the feedback. 

The number of photons received from the filament 
detectable in Fig. 13 (bottom) is ~ 21000 for an inte- 
gration time of 200 ks. This is, in principle, 10 times 
more than strictly necessary for measuring the temper- 



ature of the feature. However, the filament photons are 
spread over the entire detector and severely contaminated 
by the point source population (resolved or not). With the 
present analysis techniques, it is not possible to extract a 
workable spectrum. 



5.4. Technical remark on filtering techniques 

Multi-resolution filtering techniques give very good results 
for enhancing the signal from point source and cluster type 
emission as well as assessing the significance of the de- 
tected objects (cf Fig. 10 & 11). However, our preliminary 
processing shows that gaussian filtering is better adapted 
for enhancing very diffuse structures having a size com- 
parable to that of the detector. In a forthcoming paper 
(Valtchanov et al 2000), we shall investigate further source 
detection and extraction techniques and the possibility of 
obtaining spectroscopic information from diffuse filamen- 
tary regions. 



6. Conclusion. How can the cosmology be con- 
strained? 



We have presented detailed predictions as to the visibility 
of the cosmic structure by XMM using high resolution hy- 
drodynamical simulations. Special attention has been paid 
to physical processes — in addition to pure gravitational 
infall and shocks — which can affect the X-ray emissivity. 
The discussion was focussed on the lowest density objects, 
i.e., small groups of galaxies and cosmic filaments. 

One of the main conclusion is that the entire (or much 
of) the low mass group population will be visible with 
modest exposure times out to a redshift of one half, de- 
pending on the adopted cosmology and details of the 
physics. More specifically, for open universes, we have in- 
vestigated two effects (cooling and feedback) using sim- 
plified prescriptions. We now place them in a more gen- 
eral context: (i) Feedback from galaxies, such as super- 
nova driven galactic winds, is expected to be particu- 
larly efficient prior to group formation since the rather 
small masses of groups makes them more sensitive to 
the entropy of the gas from which they form. Energy 
and matter input to the ICM occurs throughout the life- 
time of the group, but probably but with less impact on 
their luminosities than with the preheating adopted here 
(Metzler & Evrard 1994); (ii) Further, because of their 



lower temperatures, groups are to be also strongly affected 
by cooling. This process is essentially multi-phase and be- 
cause of the condensation of cooling clouds, the remaining 
hot gas has an X-ray luminosity strongly reduced com- 
pared to what would be expected from an equivalent adi- 
abatic group; the X-ray mass fraction is also significantly 
influenced by this multi-phase evolution. As we have em- 
phasized in this work, neither of these effects are currently 
well understood or modelled. 
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The XMM group sample should then provide precious 
insights into these various mechanisms. However, a com- 
prehensive analysis of the group population should neces- 
sarily also involve a combined optical study, since galaxy 
distributions and velocities are also key parameters in un- 
derstanding the histor y of groups (an e xample is the so- 
called j3 problem, e.g. |Gerbal et al 1994 ). 

While less numerous and significantly less affected by 
cooling or feedback, large clusters — the deepest poten- 
tial wells in the universe — play an important role in our 
current cosmological constraints. For example, their num- 
ber density provides the best estimate of the amplitude of 
the primordial density fluctuations at scales around 8/i _1 



Mpc (e.g. Viana & Liddle 1996). Similarly, their evolution 
in redshift has been suggested by many authors as a way 
to measure the matter density of the universe. However, 
a convincing measurement of this has proven difficult due 
to the uncertain relation between a cluster's mass and its 
X-ray luminosity. In fact, this has driven many workers to 
use the temperature, which is more difficult to measure, in 
addition to the luminosity since the former is more closely 



related to the mass (e.g. Oukbir & Blanchard 1997). Since 
much of the uncertainty in understanding cluster luminosi- 
ties stems from feedback, the proposed observations for 
the low mass sample (along with improved modeling) will 
be able to cast light on such things as the Lx —T relation, 
and so provide much tighter cosmological constraints. 

Even more interesting, however, is the idea of using the 
spatial distribution of groups and filaments themselves to 
constrain cosmology. Even a cursory glance at Fig. 8 shows 
that the width, length and distribution of filaments is sen- 
sitive to cosmology (in fact, the dependence at z = is 
almost entirely on the amplitude of the power spectrum) . 
The difficulty is extracting this information from the more 
realistic observations shown in later figures. One possibil- 
ity is to use the location of small groups to trace out the 
filaments and apply topological measurements, such as the 
genus statistic. This has been done in the case of galaxies. 
A difficulty with this and other approaches is simply the 
weak signal of the filaments. In this paper we have shown 
that they are detectable with XMM, but confusion will 
be an issue. A promising approach to increase the signal- 
to-noise ratio is to combine multiple measurement tech- 
niques by cross-correlating galaxy surveys, X-ray maps 
and perhaps even microwave maps, to use the expected 
Sunyaev-Zel'dovich signal. This will help to eliminate the 
confusion caused by unrelated (and hence uncorrelated) 
sources. In the optical band, weak lensing analysis tech- 
niques are rapidly improving and it may soon become pos- 
sible to measure the direct gravitational signal from the 
mass associated with filaments (preferentially located be- 
tween two clusters). Having gained a better understand- 
ing of feedback and cooling mechanisms from the study 
of groups, we may then be in a position to interpret the 
observed X-ray emissivity in terms of density, tempera- 
ture and abundances and hence infer the baryon fraction 



in filaments. 
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